Critical Level Statistics and Anomalously Localized States at the Anderson Transition 
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We study the level-spacing distribution function P(s) at the Anderson transition by paying atten- 
tion to anomalously localized states (ALS) which contribute to statistical properties at the critical 
point. It is found that the distribution P(s) for level pairs of ALS coincides with that for pairs 
of typical multifractal states. This implies that ALS do not affect the shape of the critical level- 
spacing distribution function. We also show that the insensitivity of P(s) to ALS is a consequence 
of multifractality in tail structures of ALS. 

PACS numbers: 71.30.+h, 73.20.Fz, 64.60.Ak, 71.7D.Ej 
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I. INTRODUCTION 

It is well recognized that a statistical description of 
spectral correlations provides a powerful tool to study 
the disorder induced metal-insulator transition, namely, 
the Anderson transitioniii2i2ii Among several measures 
representing spectral correlations, the nearest-neighbor 
level-spacing distribution function P(s) has been exten- 
sively studied so far, where s is the energy spacing be- 
tween adjacent levels normalized by the mean level spac- 
ing A. The functional form of P(s) is closely related 
to the localization nature of corresponding wavefunc- 
tions. According to the random matrix theory^ the level- 
spacing distribution function in the metallic phase is ap- 
proximated by the Wigner-Dyson distribution, namely, 



Pw(s) = aps exp(-c (3 s" : ), 



(1) 



where (3 (= 1, 2, and 4 for the orthogonal, the uni- 
tary, and the symplectic ensembles, respectively) char- 
acterizes the universality classes. The coefficients ap 
and c/3 are determined by the normalization conditions 
J Pw(s)ds = J sPw(s)ds — 1, as a\ = n/2, a 2 = 32/n 2 , 
a 4 = 2 18 /3 6 tt 3 , ci = tt/4, c 2 = 4/tt, and c 4 = 64/9tt. 
For s <C 1, the distribution function is proportional to 
s 13 , which implies that adjacent energy levels cannot ap- 
proach each other indefinitely because of mixing between 
two extended states. Since P\y( s ) ^ 5 4 f° r sCl for the 
symplectic ensemble, the level repulsion in this ensemble 
is strongest amongst the three classes. In the insulating 
regime, mixing of quantum states belonging to adjacent 
levels can be ignored and the energy levels are uncorre- 
cted. Consequently, the level-spacing distribution func- 
tion obeys the Poissonian, 



Pp{s) = exp(-s) 



(2) 



The Poisson distribution for the insulating phase and the 
three kinds of Wigner-Dyson distribution functions for 
the metallic phase have been confirmed numerically in 
two and three dimensions 

In addition to these two types of distribution functions, 
there exists the third distribution at the critical point of 
the Anderson transitionASuLifliiiiiSi A critical wavefunc- 
tion percolates over the whole system in a multifractal 



manner reflecting no characteristic length scales, and the 
spatial profile of the wavefunction is similar to neither ex- 
tended nor localized statesii^iiii^ Since the level statis- 
tics are governed by wavefunction profiles^ the critical 
level statistics are expected to be different from both P\v 
and Pp. In fact, it has been numerically elucidated that 
the critical level-spacing distribution cannot be expressed 
by either Eq. P or @ and is scale independent and uni- 
versal. The analytical form of the critical distribution 
P(s) is of great interest. 

There exist attempts to solve this problem in a frame- 
work of the random matrix theory. In the random matrix 
theory, a matrix ensemble is characterized by the joint 
probability distribution of eigenvalues {£{}■ This prob- 
ability distribution can be regarded as the Gibbs func- 
tion of one-dimensional interacting fictitious particles at 
the inverse temperature /3(= 1, 2, 4), whose positions are 
given by Si . If these particles interact each other via a log- 
arithmic repulsion, namely, V5 n t(|£i — £j ) = — In \ si — Sj \ 
and are in the confinement potential V po t(si) = ef, the 
Gibbs function gives the Wigner-Dyson level spacing dis- 
tribution functions Eq. QJ. In order to obtain the ana- 
lytical form of P(s) at criticality, two types of ingenu- 
ities have been proposed. One is to set the power-law 
interaction Vi n t(|£i — £j\) oc \si — £j\~ 7 in the Gibbs func- 
tion, where 7 = 1 — 1 / dv and v is the localization length 
exponent>iii& This leads the critical distribution func- 
tion proportional to exp(— c^s 1+1 / dl/ ) for s ^> 1. The 
whole profile of P(s) is then conjectured to be 



p(«) 



- n'J* 



exp(- 



J f} s 1+1 l dv ) 



(3) 



However, most of numerical results presented so far devi- 
ate from this analytical formi 19 i 20 i 21 i 22 i 23 Although P(s) 
for s <C 1 is proportional to s 13 , P(s) obtained by nu- 
merical studies behaves as exp(— s/sq) with sq < 1 in 
the large s limit. Another attempt is to construct a 
random-matrix ensemble with a log-squared potential 
^pot(^i) = (ln£j) 2 /2ai 24 i 25 Using such a non-Gaussian 
random-matrix ensemble, Nishigaki has obtained a dif- 
ferential equation for the critical distribution P(s) > 2 4 Al- 
though the numerical solution of this differential equation 
well reproduces the critical distribution P(s) calculated 
by exact diagonalizations, the theory contains an ambigu- 
ous parameter a and the closed form of P(s) remains to 
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be unknown. Thus the analytical form of P(s) at the 
Anderson transition is still unclear. 

These arguments on P(s) are based on the assumption 
that all of critical states are multifractal characterized by 
an entire spectrum comprising infinitely many exponents. 
Recently, it has, however, been shown2& that anoma- 
lously localized states (ALS) 2 7i 28 i 2fl i 3n i 31 i 32 i 33 in which 
most of amplitudes concentrate on a narrow spatial re- 
gion exist with a finite probability in infinite systems 
at the Anderson transition point. ALS are brought by 
statistical fluctuations in disorder realizations and show 
no multifractality because of the characteristic length of 
large amplitude regions. It should be noted that typi- 
cal states are kept to be multifractal in critical systems. 
When studying statistical properties of physical quanti- 
ties defined at the critical point, one should take notice 
of the influence of ALS. 34 In particular, ALS make a sig- 
nificant contribution to distribution functions of critical 
quantities. If the level-spacing distribution P(s) calcu- 
lated at the critical point is strongly affected by ALS, 
there is a possibility that P(s) for level pairs of typi- 
cal wavefunctions might be totally different from that 
obtained numerically so far. It follows that the contribu- 
tion of ALS should be eliminated from a simulation result 
for checking numerically the validity of an analytical ex- 
pression of P(s) at criticality. This is because previous 
theoretical arguments are essentially based on the ran- 
dom matrix theory which does not allow the existence of 
ALS. 

In this paper, we study numerically the influence of 
ALS to the level-spacing distribution function at the An- 
derson transition. Energy levels and their corresponding 
wavefunctions at the critical point are prepared in two- 
dimensional symplectic systems described by the SU(2) 
models By employing our recently proposed definition 
of ALS at criticality, we distinguish ALS from typical 
multifractal states (MS) and construct ALS (MS) level- 
pair ensembles which contains only pairs of energy levels 
corresponding to ALS (MS). It is found that the critical 
level-spacing distribution function for an ALS level-pairs 
ensemble coincides with that for MS level-pair ensembles 
and for the entire ensemble. All characteristics in level- 
pair ensembles are absorbed by their mean level spacings. 
This implies that the elimination of ALS does not affect 
the profile of the level-spacing distribution function. We 
also show that this remarkable feature of P(s) is a conse- 
quence of multifractality in tail structures of ALS. This 
paper is organized as follows. In Sec. [HI we give a quanti- 
tative definition of ALS at the critical point based on the 
idea that ALS do not show multifractality. Level-pair en- 
sembles constructed for extracting contributions of ALS 
are also defined in this section. Furthermore, a brief ex- 
planation of the SU(2) model employed in this paper is 
given here. Calculated results are shown in Sec. IIIII An 
interpretation of our results and a condition that the crit- 
ical level-spacing distribution function should satisfy are 
presented in Sec. IIVI together with our conclusions. 



II. DEFINITION OF ANOMALOUSLY 
LOCALIZED STATES AND THE MODEL 

In order to study the influence of ALS to the critical 
level-spacing distribution function, it is necessary to dis- 
tinguish ALS from a set of critical wavefunctions. To 
this end, we employ an expediential definition of ALS 
proposed in Ref. 18, which is based on the fact that a 
typical critical wavefunction is multifractal while ALS 
are not. Multifractal distributions of wavefunction am- 
plitudes can be characterized by the box-measure corre- 
lation function G q (l, L,r)^ 

6 " T b b r 

where ^ m = £ ieb(/) l^| 2 and Hb t (l) = J2ieb r (i) l^l 2 are 
the box measures of squared amplitudes of a wavefunc- 
tion \ipi\ 2 in a box b(l) of size I and in a box b r (l) of size I 
fixed distance r — I away from the box b(l), respectively, 
Nb (or Nb r ) is the number of boxes b(l) [or b r (l)], q is 
the moment, and the summations are taken over all such 
boxes in a system of size L. We concentrate on the I de- 
pendence of G2(l, L, r) for r = I and the r dependence of 
G 2 {l,L,r) for I = 1. Denoting Q(l) = G 2 {l,L,r = I) and 
R(r) = G%{1 = 1,L, r), these two functions obey, if the 
wavefunction is multifractal, the following power lawsj2& 

Q{l)o,l d+T ^\ (5) 

and 

R(r) K r -[*f2r(2)-r(4)] j (g) 

where d and r(g) are the spatial dimension and the mass 
exponent for the qth moment of measures, respectively*^ 
We regard wavefunctions not satisfying Eqs. J5J) and 
© as ALS at the critical point. To make this defini- 
tion quantitative, we introduce variances Var(log 10 Q) 
and Var(log 10 R) from functions log 10 Q(l) = [d + 
r(4)]log 10 Z + c Q and log 10 R(r) = -\d + 2r(2) - 
r(4)] log 10 r + cr obtained by the least-square fit for a 
specific wavefunction. Non-multifractality of the wave- 
function is quantified by 

F = AVar(log 10 Q) + Var(log 10 R), (7) 

where A is a factor to compensate the difference between 
average values of Var(log 10 Q) and Var(log 10 R). The 
quantity T is small (large) when a given wavefunction 
is close to (far from) the typical multifractal state. We 
call r the atypicality of the wavefunction. Thus one can 
expedientially regard the wavefunctions with T > r^ Lg 
and F < r^ g as ALS and MS, respectively, where r^ LS 
and rj^g are criterial values of T for ALS and MS. 

Our aim is to calculate level-spacing distribution func- 
tions for ALS and MS level-pair ensembles. The ALS 
level-pair ensemble Sals (the MS level-pair ensemble 
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5ms) is defined as a set of pairs of adjacent levels 
whose corresponding wavefunctions are both ALS (MS), 
namely, 

Wri LS ) = {(£,£') : r > ri LS and r' > ri LS }, (8) 

and 

Sus(T* MS ) = {(£, e'):T< T* MS and T' < T* MS }, (9) 

where e and s' are adjacent energy levels in the critical 
region and T and T' are the atypicalities of wavefunctions 
belonging to e and e', respectively. In addition to Sals 
and SmSj the ensemble of the whole pairs of adjacent 
levels in the critical region is denoted by So, which is 
called the entire ensemble or the original ensemble. 

In the present paper, level-spacing distribution func- 
tions are calculated for critical states in two-dimensional 
electron systems with strong spin-orbit interactions 
(symplectic systems) because of the advantage of system 
sizes. Among several models belonging to this universal- 
ity class, we adopt the SU(2) model^ because scaling 
collections are known to be negligible due to a very short 
spin-relaxation length. The Hamiltonian is compactly 
written in a quaternion form as 

H = Y, act* RiAcj > (10) 

i i,j 

where cj (c t ) is the creation (annihilation) operator act- 
ing on a quaternion state vector— and ei represents the 
on-site random potential distributed uniformly in the in- 
terval [— W/2, W/2]. (Quaternion-real quantities are de- 
noted by bold symbols.) The strength of the hopping V 
is taken to be the unit of energy. The quaternion-real 
hopping matrix element i?y between the sites i and j is 
given by 

Rij = cos otij cos /3ij t° + sin 7^ sin t 1 

— cos jij sin /3ij t 2 + sin a^- cos r 3 , (11) 

for the nearest neighbor sites i and j, and = for 
otherwise. T M (/i = 0,1,2,3) is the primitive elements 
of the quaternion algebra. Random quantities Oiij and 
jij are distributed uniformly in the range of [0, 2n), and 
(3ij is distributed according to the probability density 
P(f3)d8 = sin(2/3)d/3 for < 8 < w/2. The critical 
disorder W c of this model is known to be 5.952 for the 
energy E = 1.0— 

III. RESULTS 

We calculate the level-spacing distribution functions 
P(s) for ALS and MS level-pair ensembles at the metal- 
insulator transition point of the two-dimensional SU(2) 
model. Periodic boundary conditions are imposed in the 
both directions of systems of size L = 36. The number 
of disorder realizations is 2 x 10 5 . We extract three suc- 
cessive eigenenergies from one sample with W = 5.952, 




FIG. 1: Critical level-spacing distribution function for the 
entire level-pair ensemble So (circles). Solid and dashed lines 
represent the Wigner-Dyson distribution function for the sym- 
plectic class (0 = 4) and the critical distribution predicted by 
Aronov et al.j— namely, Eq. (JSJ with v = 0.75. The inset 
shows the same plots in a semilogarithmic scale. Dotted line 
showing an exponential decay is a guide to the eye. 



which are closest to the critical energy E = 1.0. The to- 
tal number of level spacings is then 4 x 10 5 . All of these 
level spacings construct the entire level-pair ensemble Sq. 
In contrast to conventional studies of the level statistics, 
eigenstates are also required for distinguishing ALS and 
MS from the set of critical states. The forced oscillator 
methods* 3 ^ for quaternion-real matrices has been em- 
ployed to calculate unfolded eigenvalues and the corre- 
sponding eigenvectors. 

Figure shows the level-spacing distribution P(s) for 
the entire level-pair ensemble Sq. The function P(s) is 
proportional to s 4 for s C 1, while it is slightly shifted 
toward smaller s as compared with the Wigner-Dyson 
distribution. In the limit of s 1, P(s) seems to de- 
cay exponentially as seen in the inset of Fig. ^ These 
features of P(s) are consistent with numerical results re- 
ported previouslyfi^iSSiSi In fact, when we try to fit our 
data of P(s) to the critical distribution function predicted 
by Aronov et al.^ the exponent v in Eq. © is found to 
be 0.75. This value is very close to that obtained pre- 
viously by a similar analysis for the Ando models as 
an alternative model of the two-dimensional symplectic 
system,- 1 though it largely deviates from the value be- 
lieved to be correct for this universality class (y = 2.8). 
This implies that previous numerical results with large 
energy windows correctly represent the level-spacing dis- 
tribution at criticality, which differs from Eq. at large 
s values but behaves exponentially. 

We have prepared several level-pair ensembles Sals 
and Sms by choosing values of r ALS and T^s t° calculate 
the level-spacing distributions for these ensembles. The 
parameter A in the definition of T [see Eq. Q] is set to 
be 3.0, because the average of Var(log 10 R) is about three 
times larger than that of Var(log 10 Q), while the choice 
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TABLE I: List of the level-pair ensembles employed in the present work. The ensemble So is the entire (original) ensemble. 
The quantity F* represents r ALS for ALS ensembles and T^js for MS ones. A is the mean level spacing of each ensemble, and 
N is the number of level pairs included in the ensemble. 









ALS Ensemble 






MS Ensemble 




Ensemble 

r* 

A (xl(T 4 ) 
JV 


So 

8.12 
400 000 


•S'als 
0.03 
7.68 
39 100 


•Sals 
0.025 
7.74 
59154 


•Sals 
0.01 
7.99 
226 802 


•Sms 
0.01 
8.39 
32 881 


•JMS 

0.015 
8.33 
82172 


c3 
>->MS 

0.03 
8.22 
227421 
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FIG. 2: Critical level-spacing distribution functions for ALS 
and MS level-pair ensembles. Symbols A, — , x, + , and o rep- 
resent P(s) for level-pair ensembles 5*als> "Salsi <Sms, >Sms> 
and So, respectively. Filled circles show P(s) for the mixed 
level-pair ensemble Smix with F* = 0.01. The inset exhibits 
the mean level spacings for ALS level-pair ensembles (filled 
circles) and MS ensembles (open circles) rescaled by the mean 
spacing Ao for So. The abscissa indicates Tals f° r filled cir- 
cles and rjjs for open circles. 



of A does not affect our conclusions. Statistical data of 
seven ensembles used in the present work is listed in Table 

D 

Figure |2 shows the level-spacing distribution functions 
for the ensembles of S ALS , S'als > Sj^g, S^ s , and So by dif- 
ferent symbols. We should remark that the distribution 
functions are normalized and the spacing s is rescaled by 
the mean level spacing A for each ensemble listed in Ta- 
bic HI The distribut ion function for So is the same with 
that shown in Fig. 2] It is surprising that all distribu- 
tions collapse onto a single universal curve which is iden- 
tical to the conventional critical level-spacing distribution 
function for the two-dimensional symplectic class. This 
seems to conflict with the relation between the spectral 
correlation and wavefunction profiles. Naively, the func- 
tion P(s) for an ALS level-pair ensemble S'als would be 
more Poisson-like due to the localized nature of wave- 
functions. Our numerical result indicates that ALS at 
the critical point do not affect the critical level-spacing 
distribution at all in spite of the fact that the probabil- 



ity to find ALS at criticality is finite even in an infinite 
system. This implies that P(s)'s obtained by previous 
numerical works without paying attention to ALS rep- 
resent the spectral correlation for typical critical states 
and can be used for judging the validity of theories for 
the critical level-spacing distribution based on the ran- 
dom matrix theory. In the next section, we will discuss 
the reason of this remarkable property of P(s) at criti- 
cality. The inset of Fig. [3 shows mean level spacings for 
ALS and MS level-pair ensembles rescaled by the mean 
level spacing Ao for the entire level-pair ensemble So- 
The criteria T* of level-pair ensembles means r^ Lg for 
ALS level-pair ensembles and TJ^g for MS ensembles. In 
the case of ALS level-pair ensembles (filled circles), the 
mean level spacing decreases from Ao with increasing 
r ALS , while for MS level-pair ensembles (open circles) it 
increases as decreasing T^g . This implies that ALS level 
pairs mainly contribute to the level-spacing distribution 
in a small s region as demonstrated by Table [I] which is 
consistent with the fact that the repulsive force between 
two adjacent levels corresponding to localized states is 
weak. 



IV. DISCUSSIONS AND CONCLUSIONS 

Let us consider the reason why the functional form of 
the critical level-spacing distribution P(s) is not affected 
by ALS as shown in the previous section. To answer this 
question, it is important to examine carefully the spatial 
profile of ALS. Figure [21 shows squared amplitude distri- 
butions of an ALS (r = 0.100) at criticality, a typical 
MS (r = 0.009), and a localized state in the insulating 
phase (W = 15.0 and E = 1.0) of the SU(2) model of 
size L = 120. We see from the upper row of Fig. [3] that 
amplitudes of the ALS wavefunction [Fig. EJa)] concen- 
trate on a narrow region in the system, which resembles 
the usual localized state [Fig. G^c)] in appearance and is 
in contrast to the MS wavefunction [Fig.^b)]. There ex- 
ists, however, a crucial difference between profiles of ALS 
and truly localized states. The difference can be found in 
tails of wavefunctions. The lower row of Fig.[3]shows the 
gray-scale plots of the logarithm of squared amplitudes 
corresponding to the right above wavefunctions. For the 
truly localized state [Fig.JSJf)], amplitudes in logarithmic 
scale decrease with getting away from the localization 
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120 



FIG. 3: Spatial profiles of an ALS [(a) and (d)], a MS [(b) and (e)], and a truly localized [(c) and (f)] wavefunctions. These 
wavefunctions are computed for the SU(2) model of size L = 120. The truly localized wavefunction [(c) and (f)] is an eigenstate 
of the system in the insulating phase (W — 15.0 and E = 1.0). Figures (d)-(f) represent the logarithm of squared amplitudes 
by gray-scale plots. 



center (the center of the gray-scale plot), which means 
that the wavefunction decays exponentially. On the con- 
trary, the ALS wavefunction does not possess such an 
exponential tail as demonstrated clearly in Fig.^d). It 
should be noted that the gray-scale range of Fig. |3ff) 
(1CT 36 to 10- 1 ) is much wider than that of Fig. EJd) 
(10 -10 to 10~ 2 ). A remarkable feature of the ALS wave- 
function is that the amplitude distribution away from 
the peak position of the ALS [the center of Fig. Old)] is 
quite similar to that shown in Fig. |3fc). This implies 
that amplitudes in the tail region of ALS distribute in a 
multifractal manner. 

To confirm this perspective, we performed the mul- 
tifractal analysis for the amplitude distribution only in 
the tail region of ALS. For this purpose, we have ex- 
tracted a part of the eigenstate depicted by Fig. |3fa) 
within a 60 x 60 subsystem cut from the original sys- 
tem (120 x 120) so that the extracted eigenstate contains 
only a tail region of the whole eigenstate. (For cutting 
the subsystem, the original eigenstate has been appro- 
priately shifted by taking the periodic boundary condi- 
tions into account.) In order to analyze its multifractal 
properties, we have renormalized the extracted wavefunc- 
tion and calculated the box-measure correlation function 
G q (l,L,r) for q = 1, L = L s , and r — I, where L s is the 
size of the subsystem (L s = 60). Figure 0] shows that 
the calculated Gi(l,L s ,l) is proportional to l d+D2 with 
L>2 = 1-69 ± 0.01, where Z?2 = i~(2) is the correlation 
dimension, while the correlation function G\(l,L,l) for 



the whole wavefunction is not. This value of Z?2 is very 
close to the value reported so far for typical multifractal 
wavefunctions belonging to this universality class^& This 
implies that tail regions of ALS exhibit the same multi- 
fractality with that of typical critical wavefunctions4i^ 

The tail structure of ALS gives a clue to understanding 
the behavior of P(s) at criticality for the ALS level-pair 
ensemble Sals- Due to non- vanishing multifractal tails 
of ALS, quantum states belonging to adjacent ALS lev- 
els weakly couple each other through the overlap between 
wavefunctions in their tail region. As a consequence, we 
expect a repulsive force between these two ALS levels. 
The strength of the repulsive force is weaker than that 
for MS level pairs, because the average amplitude in the 
ALS tail region is small compared to the average of MS 
amplitudes. This explains a small mean level spacing for 
Sals as presented in Tableland the inset of Fig. [2] Since 
the spectral correlation is governed by overlaps of multi- 
fractal tails of ALS wavefunctions, it is reasonable that 
the critical level-spacing distribution for an ALS level- 
pair ensemble has the same functional form with P(s) 
for MS level-pair ensembles. Multifractality in tails of 
ALS does not depend on the value of r^ LS . Therefore, 
the critical level-spacing distributions for the entire level- 
pair ensemble So, the ALS ensemble Sals with any r^ LS , 
and the MS ensemble Sms with any T^g have the same 
profile, and all characteristics of level-spacing ensembles 
are absorbed by the mean level spacings. 

The fact that the critical level-spacing distribution is 
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FIG. 4: Box-measure correlation functions G q (l,L,r) for 
q = 1 and r = I for the whole eigenstate depicted by Fig. EUa) 
(open circles, L = 120) and the part of the wavefunction 
containing a tail region cut from the original ALS wavefunc- 
tion(filled circles, L = L a = 60). Dashed line shows the least 
square fit for Gi(l,L a ,l) for the extracted wavefunction. 



FIG. 5: Numerical confirmation of Eq. 115t . Solid line rep- 
resents the critical level-spacing distribution function calcu- 
lated for the entire level-pair ensemble So. Open circles and 
crosses are obtained from the right-hand side of Eq. 1151 with 
r* = 0.01 and 0.03, respectively. The inset shows the same 
plots in a semilogarithmic scale. 



invariant under changes of level-pair ensembles gives a 
condition for P(s). If the ALS criteria r ALS is equal to 
Tms (say T*), the entire level-pair ensemble So can be 
decomposed into three ensembles Sals, Sms an< ^ ^mix, 
where S m i x contains level pairs of ALS and MS, namely, 
S mix (F*) = {(£,£') : F > F* and V < T*}. These en- 
sembles have no overlaps each other, i.e., Sals H Sms = 
Sms n S mix = S mix n Sals = 4> (null set) and S = 
Sals U Sms U S m ; x . In addition to P(s)'s for Sals and 
Sms, w e have also confirmed that the distribution P(s) 
for S m ix has the same functional form with P(s) for So, 
Sals, or Sms (see filled circles in Fig.|2J. Thus, we have 



P(S) = PALS(S) = PMS(S) = Pmix(s), 



(12) 



for any value of T*. Here, P x (s) is the critical distribution 
for S x (x stands for the suffix "ALS", "MS", or "mix"). 
For T 



ALS 



1 MS' 



P(s) can be expressed as a sum of 
contributions from Sals, Sms, an d S m i x : 

P(s) = Pals(s) + Pms{s) + P roix (s), (13) 

where the distribution P x (s) is defined to be a function of 
the level spacing rescaled by the mean level spacing Ao 
for So (not S x ) and normalized as / P x (s)ds = N x /Nq, 
where No and N x are the numbers of level pairs included 
in the ensembles So and S x , respectively. [P x (s) is a func- 
tion of the spacing rescaled by A x for S x and normalized 
to unity.] Taking into account differences in the normal- 
ization conditions and the meanings of s, P x (s) is related 
to P x (s) as 



o x 



(14) 



where n x = N x /Nq and S x 
(jl4"|) . we obtain 



A x /A . Using Eqs. fT^f> - 



P(s) 



"ALS 
<^ALS 



P(s/6 ALS ) 



^P(s/5 MS ) + ^P( s /S mix ). (15) 

OMS Omix 



Note that quantities n x and 5 X depend on T*. The crit- 
ical level-spacing distribution function P(s) should sat- 
isfy the above condition for any value of T* m& Figure 
shows numerically the validity of this condition. Data 
fabricated by using the right-hand side of Eq. (|15|) for 
r* = 0.01 and 0.03 (symbols) agree quite well with the 
critical distribution itself (solid line). 

We should note here that Eq. (|15|) is not consistent 
with the exponential behavior of the critical level-spacing 
distribution function in the large s limit, which is numer- 
ically observed in the inset of Fig. ^ Thus, Eq. (|12l) [or 
Eq. (fTo|) ] implies a possibility that the profile of P(s) 
for s 3> 1 is not exactly exponential. For clarifying the 
asymptotic behavior of P(s) for s 3> 1 or the validity of 
Eq. I|12(l . further investigations are required. 

In conclusion, we have studied the level-spacing distri- 
bution function P(s) at the Anderson transition point of 
the two-dimensional SU(2) model belonging to the sym- 
plectic class by paying attention to ALS. Since ALS at 
criticality exist with a finite probability even in the ther- 
modynamic limit, quantities at the critical point might 
be influenced by ALS. Facts that ALS seem to have very 
different wavefunction profiles from typical MS ones and 
P(s) is governed by spatial structures of wavefunctions 
intimate that the level-spacing distribution is greatly af- 
fected by ALS. To examine the influence of ALS to P{s) 
at the critical point, we prepared numerically several 
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ALS (MS) level-pair ensembles which are constructed by 
pairs of ALS (MS) levels. Our remarkable result shows 
that level-spacing distributions for ALS level-pair ensem- 
bles coincide completely with those for MS and the en- 
tire ensembles, while the mean level spacing A for the 
ALS ensemble becomes smaller than A for the MS en- 
semble. Our findings imply that the spectral correlation 
for typical critical states can be evaluated even without 
eliminating ALS pairs from the original level-pair ensem- 
ble. Since ALS do not influence the function P(s), the 
critical level-spacing distribution function can be under- 
stood by an analytical argument based on the random 
matrix theory in which ALS are not considered. We have 
also shown that the property of P(s) insensitive to the 
existence of ALS is a consequence of mult ifr act ality in 
ALS tail structures. Furthermore, it has been pointed 
out that P(s) should satisfy the condition derived from 



the invariance of P(s) under changes of level-pair ensem- 
bles. Although these results were obtained for the two- 
dimensional SU(2) model, we believe that the significant 
property of P(s) is common in other universality classes. 
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